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Abstract 



The dependence of the statistics of energy dissipation on the Reynolds nuni- 



^ ' ber is investigated in an experimental jet flow. In a range of about one decade 

OO 
O 
O 
!>■ , 
f^ , is found to be independent on Re\, while the higher moments of dissipation 



of Re\ (from about 200 to 2000) the adimensional mean energy dissipation 



^^ \ show a power-law dependence. The scaling exponents are found to be con- 

sistent with a simple prediction based on the multifractal model for inertial 

^ \ range structure functions. This is an experimental confirmation of the con- 

rS \ nection between inertial range quantities and dissipation statistics predicted 



by the multifractal approach. 
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I. INTRODUCTION 

There is considerable experimental |jl| and numerical evidence that the longitudinal 
velocity difference structure functions 

S,{i) ^ {6u{m = C,{sY'H^/^ (^^j (1) 

are affected by intermittency corrections in the scaling exponents C,q which deviate from the 
Kolmogorov self-similar prediction Cg = g/3 ||^Q {Cq is a constant possibly depending on 
the Reynolds number, e is the turbulent kinetic energy dissipation, £ is a generic length 
scale and L is the integral length scale). The estimate of C,q is now available for a wide 
range of Reynolds numbers and different flow configurations. The analysis of experimental 
data show that structure functions display the scaling behavior (|I|) for sufficiently high 
Reynolds number, i?e, and that in such conditions the longitudinal exponents C,p become 
independent on Re P,^ . There is also experimental support for the scaling exponents to be 
i?e-independent also at small Re (when the scaling behavior in (|l]) is even not observable) 
if one make use of the so-called extended self similarity analysis (ESS) 0. However, the 
precise relation between this empirical result and the classical scaling behavior ([^), expected 
from dimensional arguments, is not well understood. We will not further discuss this point 
here. 

The basic fundamental property of fully developed turbulence is that the average energy 
dissipation (e) entering in (|l|) is asymptotically -Re-independent, when adimensionalized 
with large scale variables P,|^J^. On the other hand, it is well known that the statistics 
of the local energy dissipation, which assuming isotropy can be defined in terms of its ID 
surrogate 

.W = 15.(|^)\ (2) 



strongly depends on i?e, becoming more and more intermittent with increasing Re |]TO|,|TT 



This effect is also reflected in the tails of the pdf of velocity increments at very small 



separations [p^ -p^. A statistical description of the behavior of the local energy dissipation 
(0) can be obtained in terms of quantities defined in the inertial range, such as (|1]), within the 
multifractal approach. This approach, originally introduced as a phenomenological model 
for the inertial range statistics ||T5|,|T6[, has been extended to the prediction of dissipative 
scale statistics @JT^. In particular, this approach predicts a scaling behavior of the moments 
of (0) with Re ||18|, which has been recently measured in a simplified model of turbulence 

In the present paper we investigate the Reynolds dependence of the statistics of (H) in a 
experimental water jet. We have examined a series of data obtained from LDA measurements 
which cover about one decade of Reynolds numbers. The quality of the statistics allows us 
to compute high order structure functions with reasonable accuracy and to partially resolve 
the dissipative scales. Our main result is that the dependence of the statistics of dissipation 
on the Reynolds number is found to be consistent with the multifractal prediction obtained 
by assuming a fluctuating dissipative scale. 

The remaining of the paper is organized as follows. In Section |T| we introduce the 



theoretical models for the statistics of dissipation. In Section n|, the experimental set-up 



and the data analysis procedure. Section |IV| is devoted to the presentation of the results, 
whereas concluding remarks are given in Section 0. 

II. MULTIFRACTAL DESCRIPTION OF ENERGY DISSIPATION 



A phenomenological description of intermittency is given by the multifractal model |]15 
This model introduces a continuous set of scaling exponents h which locally relates the 
velocity fluctuations at scale i entering in (P with a large scale velocity fluctuation u': 

M^)~«'(^) . (3) 

The local exponent h is realized with a probability which scales with (i/L)^^^' where Z(h) 
is the codimension of the fractal set on which the /i-scaling holds. The scaling exponent of 
structure functions (^ are obtained by a steepest descent argument over h: 



C, = ini[qh + Z{h)]. (4) 

h 

The scaling region of (|^) is bounded from below by the Kolmogorov dissipative scale 77 
at which dissipation starts to dominate, i.e. at which the local Reynolds number is of order 
1: 

!M!Z)^1. (5) 

From (^ and (j^) one obtains that in the multifractal description of intermittency the dissi- 
pative scale is a fluctuating quantity, i.e. depends on the local scaling exponent h according 

to 

f ll'T \ ^i+TT 2_ 

r^ ~ L — ~ LRe,^ '"■' (6) 



where Rex = u'X/u is the Reynolds number based on the Taylor microscale A = J Ihvu''^ / {e) . 
Below the dissipative scale the flow can be assumed smooth and one can replace the 
derivative in @) with 



/faMX 



.= 15.^^j . (7) 

Assuming that the multifractal model can be pushed down to the dissipative scale, one can 
evaluate the statistics of @) by inserting (||) and (|) into (0). One ends with the expression 



0,10 



(e^) ~ {eY J d/x(/i)i?e^'[=^^''-^+^('^)l/('+'^) ~ {erRef'' (8) 



where again the integral has been evaluated by a steepest descent argument as 

'3ph-p + Z{hy 



Op = inf 

h 



(9) 



1 + h 

The standard inequality in the multifractal model (following from the exact result ^3 = 1), 
Z(h) > 1 — 3h 1^, implies for (^) ^1 = which is nothing but the request of finite non- 
vanishing dissipation in the limit Re —* 00. For p > 1, ^p < 0, i.e. the tail of the distribution 
of e becomes wider with the Reynolds number. 



Let us remark that the prediction (^ is based on the assumption of a fluctuating dissi- 
pative scale r] according to (^). If one assume, on the contrary, that dissipation enters in (|^ 
as an average quantity one ends up with a different predictions for the exponents 6p P,p!9 



Numerical simulations with a simplified model of turbulent cascade have shown that the 
exponents (|^) are indeed observed and the alternative prediction is ruled out |]2T]]. In the 
next Section, we will see that also our experimental data are in agreement with prediction 

(!)■ 

In the following we will consider data analysis of one component of the velocity in water 
jets at different Reynolds numbers. Because of the discretization of the acquisition, we 
are forced to replace spatial derivatives with velocity differences at small scales. The key 
quantity for our discussion is a generalization of (0) over a finite scale i: 

E(i)^15.[^J. (10) 

E{i) is a convenient definition of surrogate energy dissipation if the scale £ is sufficiently 
small. The average dissipative scale dependence on Reynolds number can be obtained from 
(0) and (|1]) as 

f] ~ LRel^^ . (11) 

In the analysis of experimental data it is convenient to normalize separations with fj and we 
will consider (|10|) at fixed t = i/fj. 

In the limit of very small separations (i.e. i* ^ 1) ([TOD recovers the finite difference 
representation of the ID-energy dissipation (|^) and thus, from (|) 

(Em - M ^ Re-^^^ (12) 

with the exponents 6p given by (|^). 

In the case of separations in the inertial range, from (|I]) we have {E[t)^) / {E[t)Y — 
£C2p-pC2_ Because i = i*f], using ( pJ]) we end with the prediction 

^ ^^^ ' -Re, '-'' (13) 



rs^ 



{E{e))p 



Let us remark that the exponents in ([l2[) and ([T3|) are expected to be not very different 
(they are both zero for non-intermittent turbulence) and thus a discrimination between the 
two prediction require good accuracy. 

III. EXPERIMENTAL SET-UP AND DATA ANALYSIS 

The experimental setup consists of a water jet in a closed circuit facility as shown in 
Figure |l]. A centrifugal pump moves water from a primary tank into a settling chamber 
which is equipped with a valve to damp oscillations due to the pump. A series of contractions 
leads to a 1.5 m long pipe, with an inner diameter of 14 cm which is followed by the final 
contraction (1:50 in area) to the jet. The jet (diameter d = 20 mm) exits into a large 
water-filled tank (height 30d, width 30d, length 60d) from which the water returns to the 
primary tank. The pipe, the contraction and large tank are made of perspex to allow 
optical access to the flow. At the jet exit, the flow is axisymmetric and has no swirl; 
preliminary measurements also confirmed that it is unaffected by any external forcing due 
to the pump. The jet has a top-hat velocity profile at the nozzle exit with a boundary layer 
shape parameter (defined as the ratio between the displacement and momentum thickness 
at the outlet) equal to 3.29 and a turbulent intensity equal to about 0.021 on the jet axis. 

Velocity measurements are performed by means of a forward-scatter Laser Doppler 
Anemometer (LDA) equipped with two Bragg cells. The fringe spacing is 3.416 fim and 
the measurement volume size is about 0.1 mm, 0.1 m,m,, 0.8 m,m, along the x, y and z axes 
respectively. The LDA data are randomly distributed in time; therefore, they are resampled 
by using a linear interpolation to obtain evenly spaced samples and also to provide unbiased 
statistics |2^ . The value of Rex is sufficiently large to expect the LDA noise not to affect the 
behavior of structure functions in the IR [^ . The LDA data on longitudinal and transverse 
structure functions and on scaling exponents were compared to data obtained by Hot Wire 



Anemometry in a similar jet; the agreement was good up to 8* order p2] 



Measurements were made at x/d ~ 40 {x is measured from the nozzle exit plane), where 



the flow fleld may be considered to be approximately self-preserving and isotropic |^ . The 



jet exit velocity Uq was selected so that the exit Reynolds number Re^ = Uod/u changes 
from about 2 x 10^ to 2 x 10^. As a consequence, Rex changes from about 200 to almost 2000 
at the measurement location. The number of collected samples is about 10^ in each run. 
As usually done in practice, we resort to Taylor hypothesis to transform time differences 
into space differences by means of the local mean velocity U. The number of independent 
samples, given by UTs/2L (where Ts is the total record duration), is about 10^. Probability 
density functions of the longitudinal velocity increments have been calculated at different 
values off*. The distributions indicate that the number of samples is adequate for achieving 
a closure of the integrand associated with the structure functions Sq{i) at least up to g = 8. 
In Table | we summarize some parameters of the experiments. 

Because the acquisition rate is constant for the different runs {6t = 5 x 10~^ s) the 
smallest resolved scale 6x = U6t varies with Rex- By increasing Rex we have both an 
increasing of 6x (as Re\) and a decreasing of the dissipative scale fj according to (^TJ). 
Despite these limitations, we are confldent that are able to resolve, at least partially, the 
dissipative range in all the reported runs. In the extreme case Rex = 1840 we are able 
to resolve down to £* ~ 10 which is just at the border of the dissipative range pH]. In 



Figure ||, the second and fourth-order structure functions obtained for Rex = 230 and for 
Rex = 1840 compensated with the scaling behavior (|l|) are plotted (Kolmogorov length 
and velocity scales are used to adimensionalise the horizontal and vertical axes). The large 
scale asymptotic values are also functions of the Reynolds number and, for the second- 
order structure functions, increase almost linearly as Rex', this result is in agreement with 
that simply derived by dimensional arguments on rms and Kolmogorov velocities, i.e. that 
u'"^ / {5u{riY) ~ Rex- Scaling exponents can be derived from structure functions of Figure |^ 
within the inertial range (approximately from £ ^ 80 to £ ~ 500, using a criterion based on 
3% difference from the maximum in the third-order structure function divided by €)- The 
results (C2 = 0.708 ± 0.020 and C4 = 1-265 ± 0.035 for Rex = 230) agree with those obtained 
numerically and experimentally by other authors [Q. 



IV. RESULTS AND COMPARISON WITH THEORETICAL PREDICTIONS 

In Figure |^ we plot the average surrogate energy dissipation {E{i)), computed from (p!OD 
and adimensionalized with large scale quantity u'^/L, as a function of Rex. A constant value 
of energy dissipation is expected from energy balance arguments 0. Thus, the observed 
independence of {E{i)), computed at fixed £* = 10, on Re\ is a confirmation that the 
dissipative scales are resolved. We can give a different estimation of the dissipation on the 
basis of the energy spectrum. The advantage in this case is that we can reduce the noise 
by interpolation of the spectrum down to the dissipative scales ||^J2^. The result for our 
data, also plotted in Figure |^, is consistent with the previous one thus confirming that the 
evaluation of the energy dissipation statistics from the present data is reliable. 

In order to verify the prediction for 6p of Section ||, the scaling laws ([I^ ) of the moments 



of E{i) with Rex must be computed. In Figure §, an example of this scaling, computed at 
i* = 10, is given for different values of p. We are confident that at least up to p = 3 the 
scaling exponents can be evaluated with an error on the fit less than 5%. 

The resulting exponents 6p are given in Figure ^. They have been computed at two 
different values of the separation £* = 10, 100. As expected, for separations moving from 
the dissipative range into the inertial range, the curve of the exponents 6{p) becomes flatter 
and approaches the prediction (P^. In Figure ^ we also plot the two predictions (|^) and 
(p!3|). These are been obtained by computing the structure function scaling exponents (q (^ 
as discussed above, at the largest Rex available. Given the exponents Q, the codimension 
Z{h) has been obtained by numerical inversion of @). Finally, (H) is used to predict the 
values of 9p. 

The agreement of our data with the two predictions (|^) and ([T3|) is remarkable. The 
value of 9p obtained from E{i) at i* = 10 are very close to prediction (|9|-p^. Some deviations 
are observed, starting from moment p ~ 2.5. These are probably due to the relatively large 
value oi i* = 10, which is at the border of the dissipative range and to the statistics. Indeed, 
according to (p!0|), the computation of 9^ corresponds roughly to the computation of the 6-th 



8 



order structure function. 

Let us remark again that the difference between predictions (^) and (|T3|) is based on 
the assumption of a dissipative scale fluctuating with the local energy dissipation. Our 
experimental data demonstrate that this is indeed the case. 

V. REMARKS AND CONCLUSIONS 

In this paper, the scaling of velocity increments at very small scales in a turbulent flow 
is investigated with special focus on Reynolds number dependence. Within the multifractal 
framework it is possible to derive a relation between the scaling of the energy dissipation 
statistics with the Reynolds number and the scaling exponents of the velocity structure 
functions at fixed Reynolds number. This relation is found to fit very well experimental 
data taken into a water jet over about a decade of Rex. The one-dimensional surrogate of 
energy dissipation is estimated from velocity differences at very small scale, at the border 
of the dissipative range. It is shown that our statistics changes from prediction (p!^ ) to 
prediction (p!3| ) as the separation i moves into the inertial range. The difference between the 
two predictions is the demonstration of the fluctuation of the dissipative scale. 

A Reynolds number dependence of the moments of the velocity increments also in the 
limit of very high values (-Re^ > 1000) has been recently reported by different authors [§. 
Due to the fact that the Reynolds number is expected to affect mainly the dissipative scales, 
the previous findings confirm experimentally the connection between small and inertial range 
scalings. This could support ESS (for example in respect to the enhanced range of scales 
involved), but simultaneously invalidate it (in respect of the independence on the Reynolds 
number). If these results are dependent on the still finite value of the Reynolds number or 
on the anisotropy of the flow field is left as a subject for future investigations. 



REFERENCES 

[1] F. Anselmet, Y. Gagne, E.J. Hopfinger and R.A. Antonia, J. Fluid Mech., 140, 63 
(1984). 

[2] A. Vincent and M. Meneguzzi, J. Fluid Meek, 225, 1 (1991). 

[3] A.S. Monin and A.M. Yaglom, Statistieal Fluid Mechanics, MIT Press, Cambridge, MA 
(1975). 

[4] U. Frisch, Turbulence, Cambridge Univ. Press, Cambridge, UK (1995). 

[5] R.A. Antonia, T. Zhou and G. Xu, Phys. Fluids, 12, 1509 (2000). 

[6] R.A. Antonia, B. Pearson and T. Zhou, Phys. Fluids, 12, 3000 (2000). 

[7] A. Arneodo et at., Europhys. Lett, 34, 411 (1996). 

[8] K.R. Sreenivasan, Phys. Fluids, 27, 1048 (1984). 

[9] K.R. Sreenivasan, Phys. Fluids, 10, 528 (1998). 

[10] C. Meneveau and K.R. Sreenivasan, Phys. Rev. Lett, 59, 1424 (1987). 

[11] C. Meneveau and K.R. Sreenivasan, J. Fluid Mech., 224, 429 (1991). 

[12] R. Benzi, L. Biferale, G. Paladin, A. Vulpiani and M. Vergassola, Phys. Rev. Lett., 67, 
2299 (1991). 

[13] K. Kailasnath, K.R. Sreenivasan and G. Stolovitzky, Phys. Rev. Lett., 68, 2766 (1992). 

[14] P. Tabehng, G. Zocchi, F. Behn, J. Maurer and H. Willaime, Phys. Rev. E, 53, 1613 
(1996). 

[15] G. Parisi and U. Frisch, in Turbulence and Predictability in Geophysical Fluid Dynamics, 
Proceed. Intern. School of Physics 'E. Fermi', 1983, eds. M. Ghil, R. Benzi and G. Parisi, 
North-Holland, Amsterdam. 



10 



[16] G. Paladin and A. Vulpiani, Phys. Rev. A, 35, 1971 (1987). 

[17] U. Frisch and M. Vergassola, Europhys. Lett., 14, 439 (1991). 

[18] M. Nelkin, Phys. Rev. A, 42, 7226 (1990). 

[19] J. Eggers and S. Grossmann, Phys. Rev. A, 45, 2360 (1992). 

[20] 

[21] G. Boffetta, A. Celani and D. Roagna, Phys. Rev. E, 61, 3234 (2000). 

[22] G.P. Romano and R.A. Antonia, J. Fluid Meek, 436, 231 (2001). 

[23] R.A. Antonia, T. Zhou and G.P. Romano, Phys. Fluids, 9, 3465 (1997). 

[24] R.A. Antonia and B. Pearson, Phys. Rev. E, 62, 8086 (2000). 

[25] R.A. Antonia and B. Pearson, Flow Turbulence and Combustion, 64, 95 (2000). 



11 



TABLES 



Outlet mean velocity 
Uo (cm/s) 


102 


196 


304 


487 


1060 


Red 


20400 


39200 


60800 


97400 


212000 


Local mean velocity 
U (cm/s) 


15.3 


30.4 


41.4 


72.1 


167 


Local rms velocity 
u' (cm/s) 


4.4 


9.2 


12.3 


23.8 


50.1 


Rex 


230 


308 


435 


926 


1840 


Local integral scale 
L (cm) 


5.5 


6.3 


9.4 


22.0 


41.1 


Local Kolmogorov scale 
r] (cm) 


0.0181 


0.0097 


0.0086 


0.0066 


0.0045 



TABLE I. Parameters of the experiments on the jet at x/d = 40 (except for the outlet mean 
velocity and the exit Reynolds number computed at x/d = 0). 
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FIGURES 
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FIG. 1. Sketch of the experimental setup. 
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FIG. 2. Structure function of order 2 and 4 for different values of Rex compensated with the 
scaUng behavior (ffl). The x-axis is adiniensionahzed with the Kolmogorov scale. 
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FIG. 3. Average energy dissipation computed from definition (|To| ) (+) at £* 
the energy spectrum (x) as function of -Re^. 
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FIG. 4. Rex dependence of {E{tY) normalized with {E{E)Y for p = 2/3 (+), p = 4/3 (x), 
p = 2 (*) and p = 8/3 (□) at £ = lOry. The Unes represent the best fit with a power law. 
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FIG. 5. Exponents Op obtained from the data of Figm'e for i = lOry (+) and (. = lOOry (x^ 



The continuous line is the intermittent prediction ^, the dashed line is the prediction (13). 
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